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On the accretion disc properties in eclipsing dwarf nova 
EM Cyg 
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Abstract In this paper we analyzed the behavior of 
the unusual dwarf nova EM Cyg using the data ob- 
tained in April-October, 2007 in Vyhorlat observatory 
(Slovak Republic) and in September, 2006 in Crimean 
Astrophysical Observatory (Ukraine). During our ob- 
servations EM Cyg has shown outbursts in every 15-40 
days. Because on the light curves of EM Cyg the par- 
tial eclipse of an accretion disc is observed we applied 
the eclipse mapping technique to reconstruct the tem- 
perature distribution in eclipsed parts of the disc. Cal- 
culations of the accretion rate in the system were made 
for the quiescent and the outburst states of activity for 
different distances. 

Keywords accretion, accretion discs - binaries: close 
- binaries: eclipsing - stars: dwarf novae - stars: indi- 
vidual: EM Cygni - cataclysmic variables. 



1 Introduction 

EM Cyg is an eclipsing dwarf nova which has relatively 
long orbital period (P or b — 6' i .98). Dwarf novae are 
evolved binary stars, where the Roche lobe-filling red 
dwarf component accretes matter onto the white dwarf 
component. Loosing angular momentum, matter forms 
an accretion disc. Such systems produce outbursts, 
caused by instabilities in the disc. 
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EM Cyg was discov e red b y Hoffmeister in 1928 
Mumford fc Krzeminskil (1969) investigated instabili- 
ties on the eclipse light curves and no ted, that they 



Stic ning et al 



reflect instabilities in accretion disc 
1 19821 ) detected rapid oscillations on light curves with 
ty pical time scale of 1 4.6 an d 16.5 seconds 



Bradv fc Herczegl (|19771 ). using long-term photo- 



graphic observations found outburst cycle to be 24-25 
days and noted the seasonal changes of the light curve 
and the outburst cycle 

Because EM Cyg shows standstills ([North et al 



2000h . it is classified as Z Cam type star. So EM Cyg is 



the only eclipsing dwarf no va among Z Cam type stars. 

Using radial velocities, iRobinsonl (|1974l ) estimated 
masses to be 0.7 M Q and 0.9 M for the white dwarf 
and the red dwarf respectively. These parameters cor- 
responded to the thermal-timesc ale mass tran s fer in 

~ (|l98lh 



Stover et al 



the system, and were doubtful 
found that the orbital inclination of the system is about 
6 3° and the ma ss relation close to results obtained by 



Robinson (1974) 



North et alT(|2000h found that the spectrum of EM 



Cyg in the range 6230-6650 A is contaminated by light 
from a K2-5V star (in addition to the K-type mass 
donor star). The K2-5V star contributes approximately 
16 % of the light from the system in that band and, if 
not taken into account, has a considerable effect upon 
radial velocity measurements of the mass donor star. 
The revised value of the mass ratio, combined with 
the orbital inclination i = 67°, leads to the masses 
of 0.99Af Q and 1.12Mq for the mass donor and white 
dwarf respectively. 

Recent more precise measurements of the third star 
light by Welsh et al. (2007) showed that masses to be 
M wd = l.OM and M rd = 0.77M o . 

There are several too different estimates of the dis- 
tance to EM Cyg. From ellipsoidal variations in the 
infrared band, Ijameson et al.l (|198lh found secondary 
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Table 1 28 heliocentric moments of 
(HJD-240000). (1) - Crimean data, (2) 



minima of EM Cyg 
: Slovakian data 



53947.44411 1 


54269.48212 2 


54328.53527 2 


53953.40727 1 


54272. 38822 2 


54331.44782 2 


53954.41905 1 


54299.44092 2 


54336.39064 2 


53960.38641 1 


54314.57070 2 


54339. 29649 2 


53961. 40854 1 


54315.44400 2 


54357.33454 2 


54240.38778 2 


54319.51621 2 


54361. 40679 2 


54242.42310 2 


54320.38851 2 


54364.31680 2 


54256.38635 2 


54321. 55141 2 


54366.35136 2 


54258.43000 2 


54322.42526 2 




54267.44340 2 


54327.36756 2 





220 240 260 280 300 320 340 360 380 400 420 

JD - 2454000 



(d=1.25m, Cassegrain system) of the Crimean Astro- 
physical Observatory (Ukraine) with the FLI CCD 
camera in R band with 15 sec exposure and using the 
Pupava telescope (d=280 mm, Newton system) with 
the MEADE DSI PRO CCD camera in R band with 30 
sec exposure in Vyhorlat observatory (Slovak Republic) 
during May-November, 2007 (Fig.l, see online database 
of I.Tameson et~aLl (Il98lh determined the possible dis- |http://var.kozmos.sk| ). The quantity of observational 



Fig. 1 Long term light curve of EM Cyg, obtained in 2007 
in Vyhorlat observatory (Humenne, Slovak Republic). 



spectral type to be K2V and th e distance to the sys 



ter n to be 320 pc . iBailevi (|198lf ) using K magnitudes 



tance to be in the range 28 5-429 pc with most pr obable 
value of about 352 pc. But lStiening et al.l (|1982r i, using 
private communication of Stover suggests 400 pc in his 
calculations. 

Taking into ac count the light of the third star, 
Welsh et al.1 (l2005h supposed the distance to be 450 - 
500 pc. 

Winter and Sion (2003 ) found 411 pa rsecs using the 
M v versus P or b relation (Warner 1995| ). but they de- 
termined only 210 parsecs when they fitted the IUE 
spectra. 



runs are 8 for Crimean data and 61 for Slovakian data. 

C CD images we re processed with C-Munipack pack- 
age (Hroch 1998| k including flatfielding and debias- 
ing. The procedure of the " artificial" mean weig hted 
star (jAndronov fc Baklanovll2004t Ikim et al.ll2005l ) was 



used to calculate of the instrumental magnitudes of 
EM Cyg. Using a secondary photometric standard 



from Misselt ( 19961 ). we calculated transformation coef- 



Latest paper of iGodon et alj (|2009i ) give the dis- 
tances ranged from 200 to 380 pc for different WD+disc 
models of UV spectra. 



ficients for converting the magnitudes from the instru- 
mental system of Pupava telescope: 

m-m Q = -0.027(±0.012) + 0.98(±0.024) • (V - R ) 
Moments of minima we determined using the method 
of " asymptotic parabola" ( Andronov fc Marsakoval20 06')[ 



Model fits for u ltraviolet data (jWinter fc Sion!l2003 ; 
Godon et al.ll2009l l give mass accretion rate estimates to 
be about 10 -10 M & /yr. This value lies below the upper 
limit (2 ■ 10~ 9 MQ/yr), determined bv lCsizmadia et al 
( 2008h from the O-C diagram analy sis. 



(Table 1). To achieve a higher accuracy with period de- 
termination, we used 30 moments of minima obtained 



by the previous investigators too ( 


Mumford & KrzeminskiB 


1969 


Mumford 1974. 19751 


1980l Beuermann & Pakull 


1984 


Csizmadia et al. 


2008 


). 



According to lVoikhanskaia et al.l (|1978I ). there are no 



circular polarization above about 0.3% in EM Cyg. 
Despite low Galactic latitude (+4.28°), the redden- 
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2 Observations 

Our photometric observations of EM Cyg were obtained 
during September, 2006 using the AZT-11 telescope 



The improved elements are 

T = 2437882.86059(33) + 0.290909199(84) d • E (1) 

From Fig.l, one can see that the typical interval be- 
tween outbursts is about 25 days. Duration of out- 
burst is approximately equal to the time interval be- 
tween them. 

Individual light curves show significant flickering 
which deforms the partial eclipse shape. 

On the individual light curves the effect of the pres- 
ence of a bright hot spot in accretion disc is presented. 
The hump before the eclipses is well visible, and the as- 
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Table 2 The system parameters of EM Cyg, used in the 
eclipse models. 



Parameter 



value 



reference 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

phase 

Fig. 2 Averaged light curves for two different states of 
the brightness of the system, the eclipse model fits and the 
residuals. 



cending branch of the eclipse light curve is longer, than 
the descending one. 

The strong flickering and variability with time-scales 
from several minutes to dozen percents of the orbital 
period are presented on the light curves. 

To investigate the typical patterns on the orbital 
light curves of EM Cyg and to diminish the flicker- 
ing distortions, all phase curves were divided onto two 
groups, depending on the activity state of the system. 
Within each group the average detrended light curves 
were calculated (Fig. 2). 

There is well pronounced secondary eclipse on a 
phase 0.5 on the light curve both in the outburst and 
the quiescent states. This eclipse occurs when the sec- 
ondary (red dwarf) star is behind the white dwarf and 
the disc. Our data shows that the amplitudes of the 
secondary minima in the quiescent and the outburst 
states are constant. Thus, an occulting media has the 
same properties in both states because the secondary 
star has the constant luminosity expected. 

The shape of quiescent state light curve is typical for 
the systems with significant ellipsoidality effect. 



3 Eclipse mapping of EM Cyg 

Light curves of the partial eclipse of accretion disc in 
EM Cyg allow us to reconstruct the temperature distri- 
bution in the eclipsed parts of accretion disc using the 
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eclipse mapping technique (|Hornelll985f). If we u se the 
syst em paramete r s, det ermined by North et al. ( 2000[ ) 
and IWelsh etal] (|2007T ) (i = 67°,M rd = 0.77, M wd = 
1.0), one can see that there is no eclipse of the white 
dwarf in the system (see Fig. 3 in lRobinsonl ( 1974 1). The 
eclipse light curves become more sensitive to the accre- 
tion disc structure fluctuations in such situation. To 
minimize the influence of accretion disc non-uniformity, 
we used the averaged light curves for the outburst and 
the quiescent activity states of EM Cyg, presented in 
the Fig. 2. 

Our model takes into account radiation of the 
eclipsed part of the accretion disc, the donor star light 
with ellipsoidality and gravitation darkening effects and 
the constant term, wich includes the radiation from the 
uneclipsed part of accretion disc and the third star 
light. So, we cannot resolve the third star light frac- 
tion in the constant radiation parameter but the eclipse 
mapping results are independent from the presence of 
any additional constant sources. 



Our eclipse mapping technique (|Halevinll2007l ) based 
on genetic algorithm fitting of the light curves with the 
model, where brightness distribution in two-dimensionall 
accretion disc is modeled by the set of points. The den- 
sity of such points is proportional to the brightness of 
the selected area of accretion disc. The advantage of 
this method is the smaller number of free parameters 
of the model than in the classical eclipse mapping tech- 
niques are used. 

To model of the donor star radiation we suggested 
that it fills the Roche lobe and use the gravity dark- 
ening parameter (3 to be 1. The reflection effect was 
neglected because it is expected to be small in a late- 
type companion. 

The maps of the eclipsing part of accretion disc for 
the quiescent and the outburst states have been cal- 
culated using the system parameters, presented in Ta- 
ble 2. To build the map of accretion disc we used the 
model with 200 radiating points. 

The fits and the residuals one can see in Fig. 2. In 
the quiescent state we see well pronounced pre-eclipsc 
hump, corresponded to the hot spot in accretion disc, 
where ballistic flow coincides with the disc. 
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In the quiescent state, the secondary star contributes 
of about 64 per cent of the light to the total flux from 
the system and the third star in R band. 

The secondary eclipse deep is about 5.3 per cent 
of the flux from the red dwarf and the width of the 
eclipse is about 15 per cent of the orbital period. Taking 
into account low inclination angle, the origin of cover- 
ing might be ballistic stream or extended above orbital 
plane gaseous structures. 

For our models we used corre ction for interstellar 
extinction according Ferniel ( 1975f) Er-i ~ 0.74E B -v 
and Aji ~ 3.0i?R_/. In the case of Eb-v — 0.05 we 
have An to be about 0.1 per kiloparsec. 



Table 3 Mass accretion rates (Mq /yeax) for different dis- 
tances. 





state 


distance, pc 


quiescent 


outburst 


200 


1.2- 10" 9 


5.6- 10" 9 


300 


2.2 • 10~ 9 


1.4- 10~ 8 


400 


3.8 • 10~ 9 


3.1 • 10~ 8 


500 


6.1 • 10" 9 


6.2- 10" 8 



Taking into account of interstellar extinction, we cal- 
culated the brightness distribution in the eclipsing part 
of accretion disc for the distances of 200, 300, 400 and 
500 parsecs. 

Using the brightness distribution in accretion disc 
and the blackbody approximation, we calculated the 
brightness temperature distributions and compare it 
with predictions of the stationary accretion disc mod- 
els. In Fig. 3 the radial brightness temperature plots for 
different distances are shown. The dashed lines are the 
temperature distribution predicted for the steady state 
solutions for different accretion rates. The thick solid 
line represents the critical temperature above which 
gas is ionized and supports the steady accretion regime 
|Warner|[l995T) . 

When we use the distance to the system to be 200 pc, 
the temperature distribution fit is close to the steady 
state model distribution with 6 = 3/4 for the quiescent 
state and the outburst state both. At the same time, 
the critical temperature curve lies above the outburst 
temperature distribution curve. 

In the case when the distance to be 500 pc the tem- 
perature distributions in the outburst and the quiescent 
states both are even steeper than the steady state mod- 
els what is uncommon. 

From the temperature distribution curve in quiescent 
state we can find the size of accretion disc to be about 
0.35a where a is the binary separation. This value 



Comparing the temperature distribution with the 
steady state models, we estimated the mass accretion 
rate for different distances to the system in the qui- 
escent and in the outburst states both (see Table 3). 
The mass accretion rate estimates were calculated by 
minimizing of the mean square deviation of observed 
temperature values from the steady state temperature 
distribution for different accretion rates. 



4 Discussion 

Our mass accretion rate estimates are close to values, 
determined by the previous investigators, if we choose 
the distance to be about 200 pc. A t the same tim e, the 
disc becomes ionized according to Warner ( 1995| ) crit- 
ical temperature curve during outburst for the models 
with distance larger than 300 pc. Distances about 500 
pc demand ionized gas state in the disc even for qui- 
escent states, what is improbable. So, using ionization 
criteria during outburst state, we must accept the dis- 
tances, ranged between 300 and 400 pc. 

Another estimate of the distance we can make using 
ellipsoidal modulation of the light curve in the quiescent 
state, where such effect is mostly pronounced. Using 
mean flux of the secondary from our model, we find that 
it has rriR ks 13.6. If the secondary star having mass to 
be 0.77Mq is a normal main sequence star, the spectra l 
class must be K OV. In that case my = 5.9 ( Allenl 1973h . 
V — R = 0.64 (|Johnsonlll966r ) and the distance must 
be about 459 pc. At the same time for the system 



slightl y larger than that determined by iGodon et al 
|2009h 0.3a from UV data. 



parameters we used, to fill its Roche lobe, the secondary 
must be slightly evolved. I f we use the sp e ctral class 
value K2V, determined by Ijameson et al.l (|198ll ) the 
distance estimate will be 400 pc. This va lue is the same 
as determined by Winter fc Sion ( 2003| ) using the M v 
versus P or i, relation. 

The distance of 400 pc in our models corresponds to 
relatively high accretion rate. Taking into account that 
the outburst state and the quiescent state have practi- 
cally equal durations we estimate of the mass transfer 
rate to be about 1.7 • 10~ 8 M Q /year. This value about 
an ord er higher than determined by ICsizmadia et al" 
( 20081 ) upper limit 2 • 10~ 9 M Q /year for conservative 
mass transfer. To explain of observed stability of the 
orbital period in the case of high mass transfer rate we 
can take into acc ount of magnetic brakin g effect. Us- 
ing relations from ICsizmadia et al.l (|20081 ) we find that 
in the case of system parameters from Table 2, the 
mass transfer rate which can be neutralized by mag- 
netic braking is about 2.5 ■ 10~ 8 MQ/year. So, mag- 
netic braking of EM Cyg system can explain of the mass 
transfer rate wich we have found for the case of 400 pc 
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distance. Our present result supports the suggestion of 
Csizmadia et al.l (J2008|) that these two effects neutralize 
each other and that is why we do not see any observable 
period variation. 



5 Conclusions 



In this paper, using large volume of minima timings wc 
specified the orbital period of the system. Our R band 
light curves in the quiescent state show significant ellip- 
soidality effect. Firstly for this star we used the eclipse 
mapping technique of the outer parts of accretion disc. 
We estimated the mass accretion rate for different dis- 
tances and showed, that the distance must be in the 
range between 300 and 400 pc. Using ellipsoidal vari- 
ability of the secondary star we estimated the distance 
to be about 400 pc. Also we showed, that the magnetic 
braking effect can explain of 1.7 ■ l0~ 8 MQ/year mass 
transfer rate, which has been found for 400 pc distance. 
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Fig. 3 Radial brightness temperature distribution in ac- 
cretion disc for the outburst and the quiescent states for 
different distances to the system. The dashed lines are theo- 
retical temperature distribution for the steady state flat disc 
and the thick solid line shows the critical temperature a bove 
which gas is in steady accretion regime (|Warnerlll995h . 
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